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Abstract: We revisit the method of kinematical endpoints for particle mass determination, 
applied to the popular SUSY decay chain q — * x° — * & ~ * Xi- We analyze the uniqueness 
of the solutions for the mass spectrum in terms of the measured endpoints in the observable 
invariant mass distributions. We provide simple analytical inversion formulas for the masses 
in terms of the measured endpoints. We show that in a sizable portion of the SUSY mass pa- 
rameter space the solutions always suffer from a two-fold ambiguity, due to the fact that the 
original relations between the masses and the endpoints are piecewise-defined functions. The 
ambiguity persists even in the ideal case of a perfect detector and infinite statistics. We delin- 
eate the corresponding dangerous regions of parameter space and identify the sets of "twin" 
mass spectra. In order to resolve the ambiguity, we propose a generalization of the endpoint 
method, from single-variable distributions to two-variable distributions. In particular, we 
study analytically the boundaries of the {m^(j ), rriju^} and {ma, raja} distributions and 
prove that their shapes are in principle sufficient to resolve the ambiguity in the mass deter- 
mination. We identify several additional independent measurements which can be obtained 
from the boundary lines of these bivariate distributions. The purely kinematical nature of 
our method makes it generally applicable to any model that exhibits a SUSY-like cascade 
decay. 
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1. Introduction 

The dark matter problem of astroparticle physics [1] greatly motivates the search for neutral, 
stable and weakly interacting massive particles (WIMPs) at colliders [2]. WIMPs are also 
rather ubiquitous in Beyond Standard Model (BSM) physics at the TeV scale. Unfortunately, 
at the hadron colliders of the current energy frontier (the Fermilab Tevatron and the Large 
Hadron Collider (LHC) at CERN), the process of direct WIMP production, tagged with a 
jet or a photon from initial state radiation, suffers from insurmountable backgrounds [3-5]. 
In contrast, the chances of a discovery are typically greatly enhanced in case of indirect pro- 
duction, where the WIMPs are produced in the decays of heavier, more strongly interacting 
particles. 

Since the WIMPs are usually stable due to some new conserved quantum number, they 
cannot be singly produced in collisions of light SM particles. The prototypical WIMP ex- 
ample is the lightest superpartner (LSP) (typically the lightest neutralino Xi) i n low-energy 
supersymmetry (SUSY) with imparity conservation [6] . The superpartners are produced in 
pairs and each one decays through a cascade decay chain down to the Xi WIMP, which does 

1 More recently, it was realized that many of the features of _R-parity conserving SUSY are also shared by 
other model frameworks, such as Universal Extra Dimensions (UED) [7-9], Warped Extra Dimensions [10-12], 
Little Higgs theory with T-parity [13,14], etc. 
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a) two-body (on-shell) scenario b) three-body (off-shell) scenario 

Figure 1: The generic decay chain considered in this paper: D — > jC — > jl^B — > jl^ljA. Particles 
A, £?, C and I? are new BSM particles, while the corresponding SM decay products consist of a jet j, 
a "near" lepton and a "far" lepton /J. a) In the two-body (on-shell) scenario, C is kinematically 
allowed to decay to B, which then decays to A. b) In the three-body (off-shell) scenario, C is 
kinematically prohibited to decay to B, and decays directly to A. 

not interact in the detector and can only manifest itself as missing energy (MET). Since each 
event contains two unobserved WIMPs (with unknown mass), measuring the masses, spins, 
etc. of the new particles is a very challenging task. In recognition of this fact, in recent 
years there has been an increased interest in developing new techniques for mass [15-58] and 
spin [59-81] measurements in such SUSY-like missing energy events. 

There are three basic types of mass determination methods, which are reviewed and 
contrasted in Ref. [56]. In this paper we concentrate on the classic method of kinematical 
endpoints [15]. Following the previous SUSY studies, for illustration of our results we shall 
use the generic decay chain D — > jC — ► jl^B — * jl^ljA shown in Fig. |]. Here A, B, C 
and D are heavy BSM particles, while the corresponding SM decay products are: a QCD 
jet j, a "near" lepton l~ and a "far" lepton Ij . This chain is quite common in SUSY, with 
the identification D = q, C = x[>, B = I and A = Xi-, where q is a squark, I is a slepton, 
and X\ (Xz) 18 the first (second) lightest neutralino. However, our analysis is not limited to 
SUSY, since the decay chain in Fig. p] is rather typical for other BSM models as well, e.g. 
UED [9]. For completeness, we shall consider the two different cases shown in Figs. |l](a) and 
|l](b), correspondingly. In Fig. |](a) m B < mc, so that the C — ► l^B decay is two-body. In 
what follows, we shall refer to this case as the "two-body" or "on-shell" scenario. On the 
other hand, in Fig. [l|(b) > mc and the decay C — > /^/J^4 is three-body, leading to a 
"three-body" or "off-shell" scenario [33]. In the two-body scenario, the goal is to determine 
all four unknown masses, mo, mc, tub, and rriA- In the three-body scenario, the goal is 
to determine the three 2 unknown masses, mrj, mc, and tua- Of course, the scenario is not 
known until the data are examined, thus an additional goal of our analysis will be to identify 
the particular scenario at hand. 

The idea of the kinematic endpoint method is very simple. Given the SM decay products 
exhibited in Fig. ||, form the invariant mass 3 of every possible combination, m;;, Triji n , rriji f , 

2 In the three-body scenario of Fig. [j](b), it may still be possible to extract the fourth mass ms from the 
data, e.g. by analyzing the shapes of the invariant mass distributions [28]. However, those approaches are 
quite challenging, since the shapes depend on a number of additional factors: the experimental resolution, the 
spins of the new particles, the shape of the underlying backgrounds, etc. In contrast, here we are concentrating 
on methods which use only kinematic endpoint information and are thus immune to those detrimental factors. 

3 We shall see below that the formulas simplify if we consider invariant masses squared instead. This 
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and rriju. Since and Ij cannot be distinguished on an event-by-event basis, one has to 
introduce an alternative definition of the jl distributions. For example, one can identify the 
two leptons by their charge and consider the samples {jl + } and which are experimen- 

tally well defined. This approach turned out to be very useful for spin studies [59,61,67,80], 
since spin effects are encoded in the difference between those two distributions. However, for 
mass determination, it is more suitable to use an rriji ordering by invariant mass: 

m jl(lo) = mln{m jln ,m jlf } , (1.1) 

rriji(M) = m &x-{™>jln> m jlf} • ( L2 ) 

Both of the newly defined quantities mji^ ) and rriji^hi) also exhibit upper kinematic endpoints 
( rn 'ji{u>) an d m j?(Ti)' corres P on dingly), which are experimentally measurable. Together with 
the measured kinematic endpoints m^ ax and m^ x of the mu and rriju distributions, this 
gives 4 measurements 

max max max max i-\ o\ 

m ll i m jll i m jl(lo)> m jl(hi)> \ 1 - 6 ) 

which are known functions 4 of only 4 unknown parameters (in a, ms, mc andm^). Therefore, 
by inverting those relations, i.e. solving the so called "inverse problem" at the LHC [82], one 
would expect to be able to determine the complete spectrum, at least as a matter of principle. 

However, this determination can be ambiguous, and several alternative solutions for the 
masses may emerge, as already recognized in, for example [26,82-89]. This may happen for 
one of the following reasons. 



1. Insufficient number of measurements. The four measurements (|l.3| ) may not all be 
independent from each other. Indeed, there are certain regions of parameter space 



(reviewed explicitly below in Sec. 2^) where one finds the following correlation [26] 



(m^r) 2 = (Km) 2 + KT) 2 • (1-4) 



In this case, the four measurements ( |1.3[ ) are clearly insufficient and one has to sup- 
plement them with an additional measurement. To this end, it has been suggested to 
consider the constrained distribution to^^h), which exhibits a useful lower kinematic 
endpoint 'irfjii^ii) [21]- The distribution mju^ e> ^^ is nothing but the usual rriju dis- 
tribution over a subset of the original events, subject to the additional dilepton mass 
constraint 



m mo ' 



U ^ ™ ^ ^max ( i c\ 

— < mu < m u . (1.5) 

In the rest frame of particle B, this cut implies the following restriction on the opening 
angle 6 between the two leptons [20] 

0>|, (1.6) 



distinction is not central to our analysis. 
4 See Section 2.1 below. 
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thus justifying the notation for Tn-jiue^)- In what follows, we shall therefore always 
supplement the original set of 4 measurements (1.3) with the additional measurement 



of m min 

01 rn jii(e>iy 



m max max max max min /-. 7 \ 

m ll i m jll i m jl(lo)i m jl(hi)T m jll(8>l )' \ l -<) 



so that in principle there is sufficient information to determine the four unknown masses. 
Then, we shall concentrate on the question whether this determination is unique or not, 
i.e. we shall be concerned only with discrete ambiguities. As discussed in more detail 
below in Sections [2| and these discrete ambiguities arise due to the very nature of the 
mathematical problem: the relations giving the endpoints in terms of the masses are 
piecewise- defined functions, i.e. their definitions depend on the values of the indepen- 
dent variables (the masses tra, tub, nic and mo). Since the masses are unknown, it 
is not clear which definition is the relevant one, and one must consider all possibilities, 
obtain each solution, and test for consistency at the very end. 

2. Experimental resolution. Ideally, the procedure just described would yield a single 
consistent solution. Indeed, this is what happens throughout a large portion of the 



parameter space. One should keep in mind that the measurements (1.7) inevitably 
come with some experimental errors, so that within those experimental uncertainties, 
two or more solutions are possible [26,27,83]. One specific example of this type is shown 
in Table |l] for the SPSla(a) mass spectrum, which was extensively studied in [26]. Even 
with 300 fb _1 of data at the LHC, the residual experimental uncertainties (due to the 
finite detector resolution, statistical and systematic errors, etc.) will still allow two 
solutions: a "true" and a "false" one, as shown in the table. However, this is not a 
true ambiguity in the sense that it arises simply due to limitations in the experimental 
precision. With time, the latter would be expected to improve and the ambiguity may 
eventually get resolved. For example, the statistical errors would be reduced with even 
more data 5 . 

3. Non-uniqueness of the inversion. Even in the ideal case of a perfect experiment, which 
would yield results for all the five measurements Jl.7|) with zero error bars, there may 
still be multiple solutions to the inversion problem. One of the main goals of this paper 
is to identify the specific circumstances when this takes place. In Section |3] we shall 
analyze the physical mass parameter space of SUSY and find a sizable portion in it where 
an exact duplication occurs, i.e. if Nature chooses a SUSY spectrum from that region, 
the measurements (|1.7[) will be consistent with two and only two SUSY mass spectra: 
the nominal one, plus a "fake". We emphasize the fact that the duplicate solutions we 
find yield mathematically identical values for all five experimental observables in fll.Tj ). 
Therefore, neither improvements in the experimental resolution, nor increased statistics 



It is also worth noting that ref. [26] conservatively assigned a rather large systematic error for the "ij^^iL) 
measurement, since the analytical shape of its edge was unknown at the time. Since then, the shape was derived 
in [32], so that by now the threshold measurement mjll^s.) should be considered on equal footing with the 
other measurements in (p..7|). 
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Variable 


SPSla(a) [26,83] 


SU1 [84] 


SU3 [84] 


AT ■ 1 

.Nominal 


True 


lalse 


True 


False 


True 


lalse 


m xl 


96.1 


96.3 


85.3 


137.0 


122.1 


118.0 


346.8 


m ln 


14o.U 


1 A Q 1 


loU.4 


254.0 


127.5 


155.0 


411.1 


A2 


1 /6.8 


177.0 


165.5 


264.0 


245.9 


219.0 


451.6 


mq L 


53 1 .1 


53 f .5 


523.5 


760.0 


743.6 


631.0 


899.9 


Region 


(1,1) 


(1,1) 


(1,2) 


(1,1) 


(1,3) 


(1,3) 


(1,1) 


m u 


77.0 


77.0 


77.1 


61 


100 


m jl(lo) 


298.3 


298.3 


299.6 


194 


322 


m Ti(ht) 


375.6 


375.6 


375.7 


600 


418 


^max 
m jU 


425.8 


425.8 


425.6 


609 


499 


m jii(o>f) 


200.6 


200.6 


205.1 


143 


148 


247 


214 



Table 1: Examples of "mass ambiguities" previously reported in [26,83,84]. The nominal values 
for the mass spectrum are given in the leftmost column in each case. The analysis of [26,83] for the 
SPSla(a) study point used all five available measurements (1.7), and included detector resolution 
effects and statistical and systematic errors. As a result, in the case of SPSla(a) there are two 
solutions: the "true" one is in the correct region (1,1) and is close to the nominal values, while the 
"false" one is in the wrong region (1,2), but nevertheless matches all of the observed invariant mass 
endpoints (1.7) within the experimental uncertainties. The study points SU1 and SU3 are taken 
from [84], where one requires a perfect match to only the four kinematic endpoints fll.3| ), ignoring any 

vain 
«(9>f )' 



experimental errors. In this case the true and false spectra predict different values of 



will be able to resolve this duplication. Several studies in the literature [26, 34, 83, 84] 
have already raised the issue of a potential ambiguity in the SUSY mass determination. 
Some representative examples from those works are shown in Table [l]. As already 
mentioned, the duplication found in [26, 83] in the case of SPSla(a) was simply due 
to the experimental uncertainties, and would be resolved in a perfect experiment. On 
the other hand, the duplication in the case of SU1 and SU3 found in Ref. [84] is exact, 
but relies on only four (namely, the set ( |1.3|) ) out of the five available measurements 
( |1.7[) . As seen from Table [l], the inclusion of the threshold m jiiJ l > iL-\ wu l i n principle 
resolve the ambiguity. In contrast, we use the full set of measurements (|1.7| ), including 
m^/(e>-)' and we still find exact duplication. In this sense, our findings, first reported 
in [86,87], are new, and extend the results of [26,34,83,84]. For example, we find that 
exact duplication occurs in the (2,3), (3,1) and (3,2) parameter space regions (according 
to the classification in Section ^), while the examples in Table |l] belong to regions (1,1), 
(1,2) and (1,3). It is also worth pointing out that the duplicated regions of parameter 
space that we find are not consistent with a typical MSUGRA-type scenario, which may 
explain why this problem has not been more broadly appreciated earlier. Numerical 
examples of a duplication similar to ours have previously been presented in [82], and 
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our analytical results below in Section || now help understand their origin. 



Having identified the problem of duplication in the measured mass spectrum, in the 
second part of the paper we present a new method for its solution. As already emphasized, the 
two-fold ambiguity in the spectrum is exact, so it cannot be resolved by simply improving the 
experimental precision on the kinematical endpoint measurements (1.7). Instead, additional 



experimental input is needed. One option is to consider a longer decay chain, which would 
yield several additional endpoint measurements. For example, the decay chains in Fig. |l]may 
begin with an even heavier particle (say, E), at the expense of a single new parameter (the 
mass of particle E) [27]. However, the presence of such a decay chain in the data is a model- 
dependent assumption and is by no means guaranteed. Alternatively, one may supplement 



( L7 ) with data from a future lepton collider [83] , but its existence is also an assumption and 
is by no means guaranteed. Therefore we do not consider these possibilities here. 

Instead, we concentrate on the question: What additional information, which is already 
present in the hadron collider data, can be used to resolve the ambiguity? It is important to 
realize that in very general terms, the kinematics of the decay in Fig. [I] is governed by some 
three-dimensional differential distribution 

d 3 T 



da d(5 d'y 



X 



where a, (3 and 7 are some suitably chosen angles specifying the particular decay configuration 
(see, e.g. [61]). Through a change of variables, these angles can be traded for three invariant 
mass combinations of the visible decay products in Fig. [j], e.g. mu,mji+,mji- [69], but other 
sets are equally possible, let us denote a generic such set by {mi, 1712,1713}- In place of Ql.S| ) 
one then has 

d 3 T 

. (1.9) 

dm\dm2dm-i, 

The distribution ( |1.9| ) is experimentally observable and is nothing but a three-dimensional 
histogram. It contains the full information about the decay in Fig. [I], including the parti- 
cle mass and spin information. The only disadvantage of (|0]) is that it cannot be easily 
visualized. 

In order to obtain a kinematic endpoint for some mass parameter, say mi, one then 
simply integrates over the other two degrees of freedom, and builds the one-dimensional 
distribution 

- — = / dm2dm 3 

ami J am\dm2dmz 

This, being a one-dimensional distribution, exhibits an upper endpoint m r [ iax . However, in 
the process of integration in ( |1.10| ) , one is losing a certain amount of the original information 
contained in ( |l.9[) . Some of this information can be recovered if we consider a two-dimensional 



dm,2dmz — . (1-10) 



(bivariate) distribution, e.g. in (1711,1712)'- 



d 2 T f d 3 T 



dmidi7i2 J dmidi7i2dm3 
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This, being a two-dimensional distribution, will exhibit not an endpoint, but a boundary 
line, which can be parameterized by a single parameter t as (mi(i),m2(i)). Finally, if we 
stick to the original three-dimensional distribution Ql.9| ), we will obtain a boundary surface, 
parameterized by two parameters, t\ and t2, as (mi(ii, £2), ^2(^1, ^2), ^3(^1, ^2))- Given that 
bivariate and trivariate distributions are more informative than the simple one-dimensional 
histograms, it is rather surprising that they have not been used more often in the previous 
analyses of SUSY mass determination. 

The second part of the paper is thus devoted to the analysis of bivariate distributions of 
the type ( 1.1 1| ) 6 . In particular, in Section |] (Section |5|) we analyze the boundaries of the 
bivariate distributions in terms of m^i^ and m ^i^ ( m |// an d wiy)- We show that the shapes 
of those distributions are very distinct and can be used to identify qualitatively the type of 
spectrum at hand, thus resolving the duplication discussed above. We also provide analytical 
formulas for the boundaries of the kinematically allowed regions, which can be used to further 
quantitatively improve on the mass determination (see also [89]). Clearly, fitting to a line 
would yield a better precision of determining the mass parameters than simply fitting to a 
point. What is more, we shall show that the bivariate distributions offer the possibility of 
several additional measurements, in addition to those in (|l.7| ). These are the locations of some 
special points on the boundary lines, for which we provide analytic expressions in terms of 
the masses mA, mB, mc an d mo- These special points are typically hidden as subtle features 
of the one-dimensional distributions but are transparent on the bivariate distributions which 
we are advertising here. 

In conclusion of this section, we summarize the main goals and results of this paper and 
point to the sections where those results can be found. 



Analytical solution of the inverse problem. In Section pT2| we provide analytical formulas 
which allow one to calculate directly the BSM mass spectrum tua, tub, mc and mo 
in terms of the experimental inputs (|l.7| ). Our formulas are completely general, for 
example, they are valid for both the on-shell scenario of Fig. |l|(a) as well as the off- 
shell scenario of Fig. |l|(b). In addition, they can be applied to all regions in parameter 
space. The availability of exact analytical expressions for the mass spectrum in terms 
of the observed kinematical endpoints makes numerical fitting (e.g. with a program like 
Fittino [90]) unnecessary. An important simplification in our approach is that we only 
need to consider four different cases, as opposed to the 11 cases usually discussed in the 
literature. 

In Section || we identify the complete SUSY mass parameter space where exact dupli- 
cation occurs, i.e. two very different mass spectra predict identical values for all five 
endpoint measurements (FT 



In Section || we analyze the shape of the bivariate distribution in terms of tti^qs an d 
m2 ji(hi)' ^ e identify the characteristic shape of the boundary lines of the distribution 



'Preliminary results of our work were reported in [85,87]. Similar ideas were discussed more recently 
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for each of our four parameter space regions. The shape not only allows to resolve 
the ambiguity discovered in Section ||, but also contains a lot of additional useful in- 
formation. For example, the shape analysis yields an additional measurement of an 
"edge" point, , and also allows to determine the endpoints and m™ x of the 

underlying rriji n and rriji f distributions. The analytic solution to the inverse problem 
(presented in Appendix |S|) takes a particularly simple form if we make use of these new 
measurements and consider the alternative set {m™ ax , fn^i m1 ji^ x , rn™ x }. 

• In Section || we perform a similar shape analysis of the bivariate distribution in terms 
of fn 2 jn versus m^. 

2. Analytical results 

In this section we present the analytical formulas which allow one to go from the mass 
spectrum to the experimentally observable endpoints (Sec. ^7l|) and vice versa (Sec. |2.2| ). 

Before we begin, we introduce some notation. Following existing studies in the literature 
[30,61,67,80], we shall redefine the original mass parameter space 

{m A ,m B ,m c ,m D } (2.1) 

in terms of an overall squared mass scale, m D , and squared mass ratios 7 

Rij = 9 > (2-2) 
mj 

where i,j £ {A, B,C, D} . Note that there are only three independent squared mass ratios 
in (^), which we shall take as the set {Rab, Rbc, Rcd}- However, in what follows we 
shall also make use of the other ratios, e.g. Rac, Rad and Rbd, whenever this will lead 
to a simplification of our formulas. Of course, the latter are related to our preferred set 
{Rab, Rbc-, Rcd} due to the transitivity property 

RijRjk = Rik ■ (2-3) 

Notice also the useful identity 

RijRkl = RilRkj- (2-4) 

We also require all of the mass parameters (|2.2| ) to be positive semidefinite. Our analysis 
assumes three additional absolute conditions on these parameters. 

Rab < 1 , Rac < 1 , Rcd < 1 ■ (2.5) 

This imposes a general mass hierarchy, 

< rriA < rnc < mp, (2.6) 



7 The practice of redefining the parameter space in terms of squared mass ratios is quite common in the 
literature. For example, our variables {Rcd , Rbc , Rab} exactly correspond to the variables {x,y,z} used 
in [61,67,80] and the parameters {Rc,Rb,Ra} used in [30]. 
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while for the mass of B the only constraint is rriA < vtlb- Depending on the mass of B, we 
can obtain either the on-shell scenario of Fig. |l|(a), in which rriA < tub < Trie, °r the off-shell 
scenario of Fig. 0(b) , in which mp < tub and possibly even rriB < mg . In summary, we shall 
use 

{m D ,R A B,RBC,RcD} (2.7) 
as our default parametrization of the 4 dimensional mass parameter space fl2.l|) . 

2.1 Forward formulas 

Here we list the well known formulas for the endpoints (|0]) in terms of the parameters Q2.7| ) 
introduced above. 

2.1.1 On-shell scenario 

In the on-shell scenario the kinematical endpoints are given by the following formulas: 

a = (mD 2 = m 2 D R CD (1 - Rbc) (1 - Rab); (2.8) 
' m 2 D (l - Rcd)0- ~ Rac), for Rod < Rac, case(l,-), 
m 2 D (l - Rbc){^ ~ RabRcd), for R B c < RabRcd, case (2, -), 



[rrij 



'jll 



m 2 D (l - Rab)0- ~ Rbd), for Rab < Rbd, case (3, -), 



(2.9) 



^m 2 D (l- ^/Rad) 2 , 



otherwise, 



case (4, — ); 



m jl(lo) 



a = I m jl(hi) 



"jln J ' 



m 



for (2 - Rab) 1 < Rbc < 1, case (-, 1), 

for R AB < Rbc < (2 - Rab)' 1 , case (-, 2), (2.10) 



in 



for < Rbc < Rab, 



case (—,3); 



m™{f\ , for (2 - Rab)- 1 < Rbc < 1, case (-, 1) 



m 



for R AB < Rbc < (2 - -Rab)" 1 , case (-, 2) 



m max\ ; f or < -R BC < -Rab, 



case (—,3); 



where 



(2-11) 



{m?C) 2 = ™1(1-Rcd)(1-Rbc), 
Kh X Y = m 2 D (l- Rcd) (1 - Rab) , 



max \ 



ml (1 " Rcd) (1 - ^ab) (2 - -Rab)" 1 ■ 



(2.12) 
(2.13) 
(2.14) 
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The physical meaning of the latter three quantities will become clear in the course of the 
discussion i: 
is given by 



discussion in Section ||. Finally, the endpoint rrij^jgyn^ introduced earlier in the Introduction, 



e = { m M0>t)) 2 = l^Mi 1 ~ R ab){^ ~ Rbc)(1 + Rcd) (2.15) 
+ 2 (1 - R AC )(1 - R CD ) - (1 - RcdW^ + Rab) 2 {1 + Rbc) 2 - 16R AC >• 



The physical meaning of the latter quantity will be revealed in Section || In (|2.8 - 2.15| ) we 
have introduced some convenient shorthand notation 



/„mai\2 i I 'max\% ( ' max \ 2 j ( ' max \ 2 ( 'rain \ 2 

a = {mil ) > b = { m jii ) > c = [ m ji(io) ) > d = J ' e = i m i"(e>f ) ) 

(2.16) 

for the kinematical endpoints of the mass squared distributions 8 . 
One can see that the formulas ( EDpTj] ) 

are piecewise-defined: they are given in terms of 



different expressions, depending on the parameter range for Rab, Rbc and Rod- This divides 
the {Rab ■, Rbc i Rc d} parameter space into several distinct regions, illustrated in Fig. |2[ 
Following [26], we label those by a pair of integers (Nju, Nji). As already indicated in eqs. (|2.9| - 
2.11| ), the first integer Nju identifies the relevant case for m^ ax , while the second integer Nji 
identifies the corresponding case for ( m ji^y m Ji^hi))- ^ ne on-shell case considered here, 
only 9 out of the 12 pairings (Nju, Nji) are physical, and they are all exhibited within the 
unit square of Fig. ^. The remaining two regions (5,4) and (6,4) seen in Fig. |2| correspond 
to the off-shell case and will be introduced below in Sec. |2.1.2j . Notice at this point that the 
formula for the ^ii(q > il) threshold is unique. 

Using (|2.S|) , ( gjj ) and Q2.ll ), it is easy to check that the relation ( |1.4| ), which can be 



equivalently rewritten in the new notation as 

b = a + d, (2.17) 

is identically satisfied in regions (3,1), (3,2) and (2,3) of Fig. |2|. Therefore, in these regions one 
would necessarily have to rely on the additional information provided by the measurement 
of the e endpoint ( 2.15| ). 



2.1.2 Off-shell scenario 

We now list the relevant formulas [33] for the off-shell scenario of Fig. 0(b), in which mg > 



8 Note that Ref. [26] uses a, b, c, d to label the same endpoints, but for the linear masses instead of the 
masses squared. 
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Figure 2: A slice through the {Rab, Rbc, Rcd} parameter space at a fixed Rcd = 0.3, illustrating 
the 11 parameter space regions (Nju, Nji) resulting from the piecewise-definition of the m"lj lx endpoint 



in eqs. Q and ( plgj ) and the (mf^^m™^) endpoints in eqs. ( |2.10D , ( |2.1l| ), fl2.20| ) and ( |2.21 
The four Nji regions used later in our analysis are color-coded as follows: Nji — 1 (green); Nj 
(magenta); Nji = 3 (cyan) and Nji — 4 (yellow). 

m c , i.e. R B c > 1: 



a = (mn 2 = m 2 D R CD (1 - ^R^f; 

, m D(! - Rcd)(1 ~ Rac), for Rcd < Rac, case (5, 

) 



b = 



TXIa 



m 2 D {l-^RXE) 2 , 



otherwise, case (6,—); 



m 



jl(lo) 



A = I m max 
Q ~ 1 m jl{hi) 



2 1 



m 2 D (l - Rcd)0- - Rac), case (-, 4); 



e = [m 



W>f) 



m 2 D (l - R C d)(1 - Rac), case (-,4); 
= X -m 2 D {\ - v r R A ^)^2R CD (l ~ VRm) 
+(1 - Rcd) (3 + v^Rac - ^I + ^ac + 6/R^^ |. 



(2.18) 
(2.19) 

(2.20) 
(2.21) 
(2.22) 



Notice the absence of the B index in these expressions, indicating that they are indeed 
independent of the mass of the heavy (off-shell) particle B. Nevertheless, the off-shell 
case can still be represented in the (Rbc,Rab) parameter plane of Fig. as the right- 
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most yellow-shaded region. Its left boundary is the line Rbc = lj beyond which particle B 
becomes on shell, while its upper boundary is the line RabRbc = Rac = 1> beyond which A 
is heavier than C and C becomes the LSP, which contradicts our original assumption ( |2.6| ). 
For consistency with the earlier notation (Nju, Nji) for the on-shell parameter space regions, 
we shall simply use Nji = 4 to label the single off-shell case for {c, d} = {(mji^) 2 > ( rn1 ji(hi)) 2 } > 
and Nju = 5,6 to label the corresponding two off-shell expressions for b = (m™ x ) 2 given in 
eq. This gives us a total of 11 allowed (Nju, Nji) combinations, which are all exhibited 



in Fig. H 

2.2 Inversion formulas 

Having presented all the "forward" formulas for the five kinematic endpoints a, b, c, d and 
e in terms of the mass parameters Rij and mo, we are now in position to tackle the inverse 
problem: deriving the inverse relations, which would give the mass spectrum m^, tub, thc 
and tod in terms of the measured endpoints a, b, c, d and e. Our goal will be to obtain the 
exact analytical inverse formulas for each of the relevant parameter space regions of Fig. |2[ 
Until now, the inverse relations have been derived for only 6 of the 11 regions, namely (1,1), 
(1,2), (1,3), (4,1), (4,2) and (4,3), and have never included the e measurement [26]. 

Before we begin, we need to make a decision about the following issue. In general, the 
system appears to be over-constrained, since we are trying to solve for four unknowns (in a, 
m B, m c an d ttid) in terms of five measurements (a, b, c, d and e). Therefore, for the purpose 
of inversion, we are allowed to drop one of the five measurements and use only the remaining 
four. Which measurement should we drop? This question actually turns out to be quite 
important for the subsequent discussion. 

The approach of Ref. [26] (which considered only the on-shell case of Fig. |](a)) was to 
eliminate e and use only a, b, c, and d. The reasoning was that the "forward" expression 
for e (2.15) appears to be too complicated to be tackled by analytic means. However, the 



problem with this approach is that it cannot be applied in the three on-shell regions (3,1), 
(3,2) and (2,3), where the three measurements a, b and d are not independent, due to the 
relation ( |2.17| ). Therefore, in order to obtain inverse relations valid over the full parameter 



space of Fig. g, we must make use of the m ^\\{Q > iL) measurement ( 2.15 ). For the same reason, 
we must also use the rri^ff^. measurement ( |2.10| ). Therefore, the choice of candidates to be 
omitted is narrowed down to three: a, b and d, i.e. precisely the problematic ones entering 
the linear dependence relation ( 2.17] ). 



Leaving aside the experimental issues of precision, at this point it should be clear that 
it is most convenient to drop the b measurement and always perform the inversion in terms 
of a, c, d and e. There are two important advantages of our approach: 

• Since we are never using the b measurement, the linear dependence ( |2.17| ) between a, 
b and d never becomes an issue, and the same four inputs a, c, d and e can be used in 
all parameter space regions (Nji, Nju). 

• More importantly, once we eliminate b from the discussion, we do not have to worry 
about the division of the parameter space into regions labelled by the integer Nju. 
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Instead, the full parameter space is now divided just into the four color-coded regions 
of Fig. |2], each of which is uniquely identified by the value of Nji and from now on will 
be labelled as IZn^' 

1. Region TZ\ (Nji = 1), defined by i~k AB ^ Rbc < 1 and shaded green in Fig. ||. 

2. Region 7^2 (Nji = 2), defined by Rab < Rbc < 2-R AB an< ^ shaded magenta in 
Fig-!- 

3. Region IZ3 (Nji = 3), defined by < Rbc < Rab and shaded cyan in Fig. ^. 

4. Region TZ^ (Nji = 4), defined by 1 < Rbc and shaded yellow in Fig. 

In what follows, we sometimes refer to regions 1Z\ , 7^-2 and IZ3 collectively as the "on-shell" 
region, and region Tt\ as the "off-shell" region, in reference to whether particle B is on-shell 
or off-shell. This distinction is in one-to-one correspondence with the distinction between the 
two-body scenario of Fig. ||(a) and the three-body scenario of Fig. [l|(b), respectively. Note 
that the region identification only depends on the two mass parameters Rab and Rbc- F° r 
comparison, the original endpoint method utilizing the b measurement, required all eleven 
regions of Fig. |2[ whose definitions depend also on Rcd, and one must check the solution for 
consistency in each region by trial and error [26]. Instead, we have now reduced the number 
of regions from eleven down to four. Furthermore, in Section |i| we shall show that the shape 
of the kinematical boundaries of the m 2 ,,,^ versus m 2 ,,, distribution reveal the exact region 

]l(hi) ]l(lo) o 

IZi in which the mass spectrum occurs, thus eliminating the need for trial-and-error inversion 
altogether. We consider this to be one of our most important results. 

With those preliminaries, we are now ready to present our inversion formulas which can 
be cast in the following form common to all regions: 

m\ = d (a> - 1) (pi - 1) ( 7i - 1) , (2.23) 

m\ = d {on - 1) (Pi - 1) 74, (2.24) 

m 2 c = G i (a i -l)p ili , (2.25) 

m 2 D = daiPiji. (2.26) 

where the subscript i = 1,2, 3,4 is used to indicate the corresponding (color-coded) region IZi 
of Fig. ^. The quantities G,, a^, Pi, and % are functions of the measured endpoints {a, c, d, e} 
and are region-dependent, just like the "forward" expressions for the endpoints in terms of 



the input masses (see Sec. 24). Before defining Gi, cti, Pi, and ji, we identify an ubiquitous 
combination of observables 

g = 2e - a (2.27) 

and use it in place of e, so that our starting point is the equivalent set of four measurements 
{a,c, d, g}. Then the quantities appearing on the right hand side of eqs. (2. 22-2^2^) are 
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defined by 



G 2 
G A 



g(2d-g)-2c(d-g) 
9 

g(2d-g) (d-c) 
g(d-c) + 2c(d- g y 

(g(2d-g)-2c(d-g))d 
gd + 2c{d-g) 



-d + g + ^/(2d-g)g, 



a 2 



«3 



«4 



a + Gi 
G\ ' 

a + G 2 
' G 2 ' 

a + G 3 
~ G 3 ' 

_ a + G 4 

~~~gT- 



01 = 



G 1 



71 



c 

G[ 



8 - d 



72 



(2.28) 
(2.29) 



c(d + G 3 



dG s 



73^—; (2.30) 



/?4 = 74 



d + Gj 
2Ga ' 



(2.31) 



A word of caution is in order regarding the off-shell scenario of Fig. ||(b), i.e. eq. ( p.31 ). In 
that case, particle B is far off-shell and its mass is not among the relevant parameters 
for the kinematic endpoints 9 , so that we only need to determine three unknowns: uia, rnc 
and mo- At the same time, we have one fewer independent inputs within our original set 
{a, c, d, g}, since eqs. ( p. 20 ) and ( 2.21 ) imply the additional relation 



■d. 



(2.32) 



For the purpose of inversion, in the off-shell case we chose to omit c and work only with 
{a, d, g}, which are the only three endpoints appearing in eq. ( |2.31| ). Finally, the appearance 
of the square root in (|2.31 ) should not be a problem, since in the off-shell scenario the ratio 



is bounded by 



1<^<2 
9 



V2. 



(2.33) 



The set of analytical inversion formulas (^23] - ^26| , ^2l - ^3l|) is the first main result of this 
paper. 



3. Duplication analysis 

Armed with the analytical results from the previous section, we are now ready to address the 
problem of duplicate solutions and the potential discrete ambiguities in the determination of 
the mass spectrum. Our procedure will be very simple and straightforward. We shall consider 
the four (color-coded) parameter space regions IZi in Fig. |2| one at a time, and in each case 
we shall ask the question: Is it possible that identically the same values of the endpoints 
{a, c, d, e} can be obtained from another type of mass spectrum belonging to a different 
parameter space region TZj, with j ^ il And if the answer is "yes", we shall then ask two 
follow-up questions: First, exactly in which parts of IZi and TZj does this duplication occur? 
Second, will the ambiguity get resolved by utilizing the additional endpoint measurement b 
at our disposal? 



Of course, ( |2.24 ) should only be used in the on-shell case. 
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Operationally we proceed as follows. First, it is important to realize that the "forward" 
analytical formulas of Section 2.1 provide a map Ti of the corresponding parameter space 
region TZi onto the space of values of the kinematic endpoints: 

jr. 

{mA,m B ,mc,m D }i ^ {a,c,d,e} , (3.1) 
or equivalently, using the reparametrization ( |2.7| ) 

{m D , Rab, Rbc, RcD}i ^ {a, c, d, e} . (3.2) 



Similarly, the inverse formulas from Section |2.2| provide an inverse map T i 1 from the space 
of kinematical endpoints back onto the mass parameter space: 

FT 1 

{a, c, d, e} i ~ ^ {m D , R ab ,Rbc, RcD} j ■ (3.3) 

The composite of the two maps ( [3.2j ) and ( |3.3| ) for i ^ j, is a transformation 

Tij = TJ 1 ■ Fi (3.4) 

relating parameter space points belonging to two different regions, TZi and TZj, yet resulting 
in identical kinematical endpoints {a, c, d, e}: 

T 

{m D ,R AB ,RBC,RcD}i^ {m' D ,R' AB ,R' BC ,R' CD }. . (3.5) 



The transformation Tij described in ( |3.5| ) will serve as the basis of our duplication analysis. 
The exact analytical formulas for this mapping can be trivially obtained from our analytical 
results above in Section but are rather lengthy and we shall not present them here explicitly. 
However, we note that in the three on-shell cases i = 1,2,3 they have the generic form 

r 'ab = Iab (Rab ; Rbc ) (3.6) 

R'bc = fBc(RAB,RBc), (3.7) 

R'cd = fcD (Rab , Rbc > Rcd ) , (3.8) 

m' D = m D f D (R A B,RBC,RcD), (3.9) 

where /ab, /bc> fcD and fr> are the functions defining the transformation Tj. One im- 
portant feature of the Tj map (|3.6| - p.9D is that it transforms the 2-dimensional subspace of 
dimensionless parameters {Rab, Rbc} into itself. Notice that Rab and Rbc are precisely 
the parameters definitng the four regions TZi in Fig. |2[ Therefore, for the purposes of our 
duplication analysis it is sufficient to consider the simpler transformation of 

{Rab, Rbc }j ^ {R'ab^'bc}, (3-10) 



given by eqs. Q3.6| ) and (fT?]) only, instead of the more general mapping fl3,5| ) given by all four 
eqs. (PH3T9|). 
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Figure 3: The maps T 13 : Tlx 
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72.3 (bottom two panels), 



7?-3 (top two panels) and T23 
which are implied by eq. (3.1C). In both cases the target region 7^.3 is shaded in cyan. Under T13, the 
green-shaded region ABD in the top left panel transforms into the green-hatched region A'B'D' of 
the top right panel. Under T23, the magenta-shaded region BCD in the bottom left panel transforms 
into the magenta-hatched region B'C'D' of the bottom right panel. In both cases, the transformed 
(primed) region falls completely within the boundaries of the intended target (7^.3). 



We are now ready to answer the main question posed at the beginning of this section: 
does a consistent mapping (330) exist for some pair of regions IZi and IZjl Note that the 



transformation ( |3.10| ) is not necessarily always well defined: consistency requires that the 
obtained values of {R'ab, R'bc) j belong to Region IZj, which is not automatically guaranteed 
and must be explicitly checked. To put this in more formal terms, we are only interested in 
those cases where the intersection of the image of region IZi under the transformation TJj and 
the intended target region IZj is a non-empty set: 

{Tn (Hi)} n {n,} £ 0. (3.11) 

In order to find all such occurrences, we consider all possible transformations Tij with i ^ j 
and enforce the consistency check ( [3,111 ), 

We begin with the on-shell case = 1,2,3), where there are 6 possible mappings 
Tij. For the purposes of finding the duplicated portion of parameter space, it is sufficient 
to consider only 3 of them, which for convenience of illustration we choose as T13, T23 and 
T21. The corresponding results are shown in Figs. [3| and |I[ Fig. ^ shows the effect of the 
transformation T13 : 1Z\ 1 — ► IZ3 (top two panels) and T23 : IZ2 i — ► ^3 (bottom two panels), 
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Figure 4: The same as Fig. [jj, but for the map T21 : TZ-2 1 — ► Hit where the intended target is the 
green-shaded region 72,i. Under T21, the magenta-shaded region BCD in the left panel transforms 
into the magenta-hatched region B'C'D' of the right panel. The image B'C'D' has no overlap with 
its intended target IZi , except along the BD = B'D' boundary, which is left invariant under the T\% 
transformation. 

while Fig. |] shows the map Ti\ '■ ~R-2 1 — ► T^i- In both figures, the color-shaded areas in the 
left (right) panels exhibit the original regions IZi (the intended target regions Rj). The cross- 
hatched areas in the right panels depict the actual image Tij{1Zi) of the TZi region under the 
transformation Tij. For example, in Fig. ||| T13 maps the whole green-shaded region ABD on 
the left into the green-hatched region A' B'D' on the right, while in Fig. ^ (Fig. |||) T23 (T21) 
maps the whole magenta-shaded region BCD on the left into the magenta-hatched region 
B'C'D' on the right. In accordance with ( |3.11| ), duplication occurs whenever the right panels 
in Figs. H and |] exhibit an overlap between the cross-hatched area of the image and the solid 
color-shaded area of the intended target. We see that duplication occurs in the case of T13 
and T23, but not for T21, although in the latter case points which are on opposite sides, but 
close to the boundary line BD will give rather similar values for the measured kinematic 
endpoints {a, c, d, e}. 

At this point one may wonder whether the result of Fig. ^ is sufficient to prove the 
existence of duplication. Indeed, Fig. || tells us nothing about the remaining two parameters 
Rcd and mo and more specifically about their transformed values R'cd and m' D under the 
mappings T13 and T23. Duplication will in fact not occur, if R' CD and m' D turn out to be 
unphysical, for example, if R' CD < 0, R' CD > 1 or {m' D ) 2 < 0. Unfortunately, a closer 
inspection of ( |3.8[ ) reveals that 

fcD (Rab , Rbc i Rcd = 0) = 0, (3.12) 
Jcd (Rab , Rbc i Rcd = 1) = 1, (3.13) 
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for any values of Rab and Rbc, so that Rcd is always consistently mapped within its 
definition region. Similarly, we find no problem with eq. (3.9). Therefore, the duplication 
examples shown in Fig. [3| truly represent a problem. 

We then perform a similar analysis involving the off-shell region IZ^ and find no occur- 
rences of duplication, which is not surprising, since the off-shell case is more restricted, due 
to (|2.32| ). Therefore, Figs. || and ||| already provide the final answer to the first question posed 
at the beginning of this Section: which portions of the mass parameter space (|2.7| ) exhibit 
exact duplication? We can summarize our result as follows: 

For every point with Rab < Rbc < 1 (i- e - i n region 1Z\ or T^) and arbitrary 
values of Rcd and tud, there exists another parameter space point with Rbc < 
Rab < 1 and certain (in general different) values of Rcd and mn, which would 
result in identical predictions for all four kinematic endpoints {a,c,d,e}. 

The reverse statement is not true: not every point with Rbc < Rab < 1 (i-e. in region 
IZ3) is subject to duplication. Referring to the right panels of Fig. ||, only the cross-hatched 
portions of the cyan-shaded region IZ3 are duplicated. 

Having found duplication examples for the limited set of measurements {a,c, d, e}, it is 
now time to ask whether the additional measurement of the b kinematic endpoint will help. 
We find that, as might have been expected, whenever the b measurement is independent of 
the others, the duplication goes away. Unfortunately, as already mentioned in the discussion 
following eq. ( 2.17| ), in the three subregions (3,1), (3,2) and (2,3), b is not an independent 



measurement, and thus the duplication will persist even for the full set of 5 measurements 
{a, b, c, d, e}\ In terms of the subregions of Fig. |, the two cases of duplication found in Fig. ||| 
can then be summarized as 

(3.1) (2,3), (3.14) 

(3.2) ^ (2,3). (3.15) 

As long as the original parameter space point and its image belong to these particular sub- 
regions, the resulting two sets of endpoints Ql.7| ) will be identical. 

We caution the reader not to get the impression from ( [3.14 ) and ( 3.15 ) that every pa- 



rameter space point in regions (3,1), (3,2) and (2,3) is duplicated with something. Recall 
that the boundaries of the dangerous subregions (3,1), (3,2) and (2,3) depend on Rcd, thus 
the range of Rcd values resulting in duplication will now be restricted. Therefore, in the 
presence of the additional b measurement, our previous statement about duplication is now 
modified as follows: 

For every point with Rab < Rbc < 1 and any mjj, there exists a range of Rcd 
for which exactly the same values of the five kinematic endpoint measurements 
{a, b, c, d, e} can also be obtained from a different parameter space point with 
Rbc < Rab < 1 and some other (generally different) values of Rcd and mo- 
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Figure 5: The minimum value R™^ (Rbc , Rab) required for duplication, as a function of Rbc and 
Rab- The white asterisks (circles) mark the duplicate pair of points P31 and P23 (-P32 and P23) m 
Table |. 



We shall now describe the duplicated parameter space implied by ( |Og ) and ( ^15|) a bit 
more quantitatively. As we just mentioned, for any given point in the (Rbc, Rab) plane, 
there may exist a range of values for Rcd which would cause duplication. Let us denote the 
minimum and maximum values of that range by Rq£) and -RJ^ff , correspondingly. Clearly, 
both R^D and Rq D x are in general functions of Rbc and Rab- Then, the "duplicated" 
parameter space can be simply described as the set of all points {Rab, Rbc, Rcd}, which 
satisfy the two inequalities 



Rcd (Rbc, Rab) < Rcd < Rcd x {Rbc, Rab)- 



(3.16) 



If, on the other hand, the values of Rbc and Rab are such that duplication does not occur 
for any value of Rcd, we can simply take RIq 1 ^ = R(j D x , resulting in a set of zero measure 
for (|3~Tcj ). 

Now, in order to delineate the duplicated parameter space, we only need to supply its 
boundaries Rq^ (Rbc , Rab) and R™jf (Rbc , Rab) ■ Our analysis reveals that within the 
duplication region we always find 



Rcd x (Rbc, Rab) 



1, 



(3.17) 



while the function R™^ (Rbc , Rab) is plotted in Fig. ||. Duplication does not occur in the 
uniformly red region in the upper left corner, so there we choose to plot Rq^} = R^c'd = 
in accordance with our convention. Within the rainbow-colored region in Fig. [5], duplication 
will exist for any value of mr>, as long as Rcd is larger than the Rq£) value shown in the 
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figure, i.e. for 

BgjOKRBC, Rab) < Rod < 1, (3.18) 



where we have made use of ( 3.17 ). Fig. [| reveals that the duplication region is typically 



characterized by a rather high 10 value of Rcd = m c /m D . This implies that in order to have 
duplication, particle D cannot be too much heavier than particle C. This situation does 
not often arise in typical SUSY models, where D is a squark q and C is the second-lightest 
(wino-like) neutralino x°>- ^ n models with a high SUSY breaking scale like SUGRA, the 
Renormalization Group Equation (RGE) running tends to split the squark and electroweak 
gaugino masses, so that the hierarchy mo ~ mc is rather unlikely. On the other hand, 
minimal UED models predict a rather degenerate spectrum, since the mass splittings arise 
mostly at the loop level, so that mp ~ mc is rather natural in this case [8]. 

As an illustration of the whole duplication discussion so far, let us now choose two specific 
examples of duplicate mass spectra, one for the case of ( |3.14| ) and another for the case of 
Q3.15 ). The corresponding input masses and mass ratios, as well as the resulting kinematic 



endpoints, are shown in Table ||. The first five kinematic endpoints shown in the Table were 
already discussed, while the rest are new and will be introduced below in Sections |4] and |5[ 
As an application of our previous results, let us outline our procedure of selecting each pair 
of study points in Table [| Let us start with the case of (|3Tl ). Since we know from Fig. |5] 



that the whole region (3,1) is duplicated, it is convenient to first choose the point from that 
region. We select nice round numbers like Rab = 0.4 and Rbc = 0.8. This choice is indicated 
in Fig. [| with the white asterisk inside region 111. Then Fig. | shows that = 0.686, 

therefore we choose a somewhat larger value: Rcd = 0.7. This choice of Rab, Rbc and Rcd 
already guarantees duplication for any value of mp, and we choose mo = 500 GeV (another 
nice round number). The resulting masses rriA, itib and mc can be readily computed in terms 
of rri£) and the mass ratios. We call the resulting spectrum "study point -P31", which is listed 
in the second column of Table ^. Given P31, one can use the transformation T13 to obtain 
the matching spectrum in region (2,3), which is listed in the third column of Table |2| under 



the name of "study point ^23"- In the case of ( 3.15 ), we follow a similar procedure, except 
we start with a point in region 71-2 (indicated with a white circle in Fig. |5|) and then use the 
T23 transformation to obtain the corresponding point in region (2,3). The two resulting mass 
spectra (called P32 and P23) are given in the fourth and fifth column of Table ^, respectively. 

In the on-shell case the parameters {Rab, Rbc, Rcd} belong to a unit cube, due to the 
restrictions (|2.5|). The volume of the unit cube is 1. Fig. || then allows us to calculate the 
volume fraction of this unit cube which corresponds to a duplicated parameter space region. 
The result that we find is 0.158. Then one might be tempted to say that if new physics like 
supersymmetry or UED with a cascade decay of the type shown in Fig. [j](a) is discovered 
at the LHC, there would be roughly a 15.8% probability that endpoint measurements alone 
would result in a duplicate spectrum, even under ideal experimental conditions. However, it 
is rather difficult to justify such probabilistic statements, since they are not invariant under 



10 The minimum value of Rcd that we find over the whole parameter region in Fig. |B] is Rcd = 0.4 and is 
found at Rbc = 1, Rab = 0. 
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418.33 


1083.10 


447.214 


554.133 


m D (GeV) 


500.00 


1172.57 


500.00 


610.443 




0.400 


0.920 


0.200 


0.487 


Rbc 


0.800 


0.777 


0.400 


0.390 


Rcd 


0.700 


0.853 


0.800 


0.824 


T>min 
^CD 


0.686 


0.845 


0.774 


0.800 


m max (Q e y) 




145 


310 


m™ x (GeV) 




257 


369 


™» ( Gev ) 




122 


149 


Km) ( GeV ) 




212 


200 


m^ >f j (GeV) 




132 


248 


m max ( Ge y) 




212 


127 


200 


183 


m$ (GeV) 




190 


112 


126 


115 


(GeV) 




122 


212 


173 


200 


m fi(eq) ( GeV ) 




NA 


122 


149 


149 


m iK(+ )(0) (GeV) 




226 


240 


214 


230 


m jii(+)( a t) (GeV) 


Vt 


263 


257 


374 


369 


m jH(+ )(a on ) (GeV) 




257 


257 


369 


369 


m iH( _)(a n) (GeV) 




190 


193 


355 


360 


m iH ( + )(a fr) (GeV) 




256 


243 


372 


367 



Table 2: Two examples of exact duplication as implied by ( 3.14 ) and ( 3.15 ). The pairs of study 
points P31 and P23, as well as P32 and P23, exhibit identical values for all five kinematic endpoints 
™i™ ax , rn™(^, to™^), ™™u x and m™^ e>1L y Point P31 belongs to Region Tlx, point P 32 belongs to 



Region TL^ , while points P23 and P' 2Z belong to Region 7^3. In Fig. S, the duplicate pair of points 
P31 and P23 (P32 and P 23 ) is marked with white asterisks (white circles). The second-to-last and last 
blocks in the table contain the endpoint measurements which are available from the two-dimensional 
distributions ( m |/(/ ), TO |/(hi)) an( l ( m U' m 'jii)> discussed below in Sections |]and[^, correspondingly. 



reparametrizations, and furthermore, they depend on the assumed (usually uniform) prior 
for the probability distribution of new physics models in mass parameter space. Any given 
model of SUSY breaking, for example, would select a preferred parameter space within the 
unit cube, and may lower or increase this naively calculated probability. The generic mass 
degeneracy in UED, on the other hand, would prefer the region Rab ~ Rbc ~ Rcd ~ 1, 
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and the duplication is much more likely. The important result from our point of view is that 
there exists a non-vanishing duplication region, and this fact alone is sufficient to motivate 
us to look for alternative methods for mass determination, which we shall undertake in the 
following two sections. 



4. Kinematic boundary lines for the Tn/jn^ versus Tn^,^ distribution 

In this section we shall analyze the shape of the two-dimensional invariant mass distribution 

d 2 r 



dm Mo) dm Mhi) 



(4.1) 



which we imagine plotted as either a scatter plot or a two-dimensional histogram with m |m ) 
on one axis and ^j^i) on 

the other 11 . The purpose of our analysis will be twofold. On 
the one hand, we shall be interested whether we can use the shape of this two-dimensional 
distribution to resolve the mass spectrum duplication problem encountered in the previous 
section. But more importantly, we shall investigate what additional kinematic endpoint 
measurements besides those already considered in ( |1.7| ), may become available in this case. 

Recall that the variables rriji^ and rriju^ were introduced in eqs. ( |lTl] ) and ( |1.2| ) as a way 
to deal with the ambiguity in the experimental identification of the "near" and "far" leptons 
1^ and Ij in Fig. ||. It is therefore not very surprising that the shape of the m-ji^ versus 
m ji(hi) distribution ( |4.1| ) that we are interested in, is very closely related to the corresponding 
m 2 , versus m 2 -, distribution 

d 2 r 



dm 2 -! dm% 



(4.2) 

In principle, both distributions ( |4.1| ) and ( |4.2j ) depend not only on the mass spectrum, but 
also on the spins and on the chiralities of the coupling constants of the particles A, B, C and 
D involved in the cascade [67,68,80]. However, the location and the shape of the boundary 



lines in the scatter plots fl4.1| ) and (4.2) are determined purely by kinematics, and do not 
depend on the spin and type of couplings. To the extent that we are only interested in these 
boundary lines, it is therefore sufficient to ignore spin effects and consider only pure phase 

are in 



space decays, in which case the analytical results for the distributions ( [4.2|) and (4.1) 
principle already available [30]. From now on, the term "shape" will therefore refer to the 
location and shape of the boundary lines, and will otherwise have nothing to do with the 
probability density of the two-dimensional distributions such as (^J) or (4.2). 



In the on-shell case, the shape of the fnji n versus Tnji invariant mass distribution (|47 
is extremely simple, and is illustrated in Fig. |6[ The scatter plot in the (m|^,m|^ ) plane 
fills the right-angle trapezoid ONPF, whose corner points are defined as follows. Point O is 



In practice, as we shall demonstrate below, it may be more convenient to plot the linear masses m^ui^ 
and rriji(hi) on the two axes, but use a quadratic power scale instead of the conventional linear or logarithmic 
scales. 
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Figure 6: The generic shape ONPF of the bivariate distribution (4.2) in the (m^ i m ji f ) plane 



simply the origin of the coordinate system. Point N (for "near") lies on the fn 2 - ln axis, and its 
coordinate is nothing but the maximum possible value of the jet-near lepton invariant mass 



n = «C) 2 = m 2 D (l- Rod) (1 - Rbc) 



(4.3) 



which was already introduced in eq. ( [2.12 ). Similarly, point F (for "far") lies on the axis, 
and its coordinate is nothing but the maximum possible value of the jet-far lepton invariant 



mass 



(4.4) 



/ = [ m Tif X ) = ™d (1 - Rcd) (1 - Rab) , 

which was already defined in eq. (2.13). Finally, the point P is the most important of the 
four corners, since it defines the actual shape of the trapezoid, once points and F are fixed. 
The coordinates of point P in the ( m ji n i m 'ji ) plane are (n,p), where n was already defined 
while p is a new quantity: 



m 



(p) 

V = [rn jlf 



fR BC = m 2 D (l- R CD ) Rbc (1 - Rab) 



(4.5) 



Since N and P share the same iTVji n coordinate n, point P always lies directly above point 
N. At the same time, the definition of p implies that 



P<f 



(4.6) 



so that point P always lies lower than point F, as illustrated in Fig. H Finally, in Regions 
7^2 and TZ^, there is one more special point, Q, which can be seen in Fig. ||: it is the point 
where the FP side of the trapezoid intersects the 45° line iTVji n = m2 ji } ■ 
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of point Q are equal by definition, and are given by 



m Me q) y = ml(l-Rcn)y^, (4.7) 



which is nothing but the quantity previously defined in eq. ( 2.14| ). The four quantities n, p, 



f and q just introduced are not all independent, but obey the relation 

- = 1 + t—?- . (4.8) 
q n 

With those conventions, the trapezoid ONPF can be equivalently defined through the para- 
metric equation of the boundary line segment FP. A convenient choice for the line parameter 
is the running value of mji n ■ Then the parametric equation of the line FP is given by 

FP: m ] lf (m] ln ) = f-t-Pm 2 jln . (4.9) 

In terms of this parametrization, the three special m|j values introduced in Fig. ^ are given 
as follows: 

/ = 74,(0), (4.10) 

p = m 2 jlf (n), (4.11) 

q = m 2 jlf (q). (4.12) 



The last equation is exactly the relation fl4.q ). 

The color-coded regions in Fig. ^ show the allowed locations of point P, and are in one- 
to-one correspondence with the color-coded parameter space regions of Fig. | (in both cases 
we use the same color coding). This correspondence is most easily seen as follows. First, note 
that the two white areas in Fig. |6| are not accessible to point P. The region with m|j > / 
is forbidden due to the relation fl4.6[ ). Similarly, the white triangular area near the origin, 
defined by 

^i f <f- m k ( 4 - 13 ) 



is also not allowed, which can be seen by using the inequality 

p>f-n, (4.14) 



following from the defining relations (fO-1.5). Therefore, point P must belong to one of 



the three colored regions in Fig. g. As can be seen from the figure, these three regions 
are distinguished based on the value of n relative to / and p: recall that eq. J4.6D already 
determines the hierarchy p < /, so that for n there are only three possible options: n can 
be smaller than p, n can be larger than /, or n can fall in between p and /. Let us consider 
each case in turn. 
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1. The case n < p < f. Point P then lies somewhere within the green-shaded area in 
Fig. ^. Using (4^-4^), it is easy to see that the conditions n < p < f imply 



n<p< f =^ - — < R BC < 1 , (4.15) 

2 - Kab 

which was precisely the defining relation for region IZi in Fig. ^. Therefore, in Fig. |6| 
we have labelled and color-coded the area with n < p < f to match the notation for 
region TZ\ used in Fig. 

2. The case p < n < f . In this case point P would belong to the magenta-shaded triangular 
area in Fig. |6|. The conditions p < n < / now imply 

p<n< f => R AB < R BC < t> — d — ' ( 4 - 16 ) 

which is the definition of region IZ2 in Fig. ^. Once again, we label and color-code this 
region to match the notation used in Fig. |2[ 

3. The case p < f < n. Now the point P should fall somewhere within the cyan-shaded 
semi- infinite rectangular strip in Fig. ||. Using ( f4.3ff4.5| ), the constraints p < f < n now 
translate into 

P< f<n => 0<R bc <R A b , (4.17) 

which is the definition of region IZ3 in Fig. ^, again justifying the notation and color- 
coding used in Fig. |6|. 



Unlike the three on-shell cases just discussed, the off-shell scenario of Fig. |l|(b) should be 
handled with care, since the "near" and "far" lepton distinctions become meaningless in that 
case. Nevertheless, the off-shell scenario can still be represented in Fig. ^, and in fact this 
representation is unique: there is a single allowed location for point P at n = f and p = 0. 
In Fig. H this unique location is indicated with a yellow-shaded circle, which corresponds to 
the whole yellow-shaded region IZ4 in Fig. ||. In other words, in the off-shell case we can 
randomly assign "near" and "far" labels to the two leptons in each event, and then the shape 
ONF of the resulting ( m 'ji n i m 'ji f ) scatter plot will be an isosceles right triangle. 

From the preceding discussion it should be clear that the two-dimensional theoretical 



distribution (4.2) contains a great deal of useful information: its shape uniquely identifies 
the on-shell parameter space region IZi, and yields the four measurements {n, f,p, q} given 
in eqs. (fO -4.5) and ((Oj) instead of the usual two {m™^ and m™ 1 ^^)- Ideally, one would 



like to preserve and subsequently extract this additional information from the experimentally 



observable two-dimensional distribution (4T) as well. We shall now show that this is in fact 
possible, using the simple intuitive understanding of the shape exhibited in Fig. |(| 

The key is to realize that the reordering (1.1) and (p..2j ) of the (m^,m^) pair into 
a ( rn 'ji(io)' m 'ji(hi)^ P a ^ r m geometrical terms simply corresponds to "folding" the trapezoid 
ONPF in Fig. || along the 45° line fn 2 ji n = m2 ji s - This procedure is shown pictorially in 
Fig. ^, where for illustration we use an example from region i.e. p < f < n. Panel 
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Figure 7: Obtaining the shape of the rn^n^ versus m^,^ bivariate distribution by folding the m 2 ^ 
versus Tn 2 lf distribution across the line m| ; = . This particular example applies to region IZ3 . 
For the other three regions, refer to Figs. |j|(a), ||(b) and ||(d). 

(a) shows the trapezoidal shape of the original fn?ji n versus invariant mass distribution 
from Fig. ||. Now suppose that we want to convert this (m^,m^) scatter plot into a 
( m jZ(io)> m2 ji{hi)) sca tter plot, simply by reinterpreting the rn 2 ln axis as m 2 ^^ and the m|> 
axis as m^u^y From that point of view, the trapezoid ONPF in Fig. ^ divides into two 
adjacent regions: OQF (blue-shaded) and ONPQ (red-shaded). Within the blue-shaded area 
OQF we have m 2 ln < rn 2 } , so that the coordinate pair {m 2 ln , m|, ) can be directly identified 
with {fn 2 i^ y m2 ji(hi))' Thus the blue-shaded area OQF in panel (a) remains unchanged and 
appears identically in panel (b), where it is marked with a blue cross-hatch. In contrast, 
within the red-shaded area ONPQ of panel (a), the coordinates rn 2 ^ and m-ji, are in the 
wrong order, and need to be reversed when going to the ( m2 ^i y m 'ji(hi)^ sca tter plot of panel 

(b) . In layman terms, this reversal corresponds to "folding" the trapezoid ONPF along the 
45° line OQ, as shown in Fig. [?]. The resulting image ON'P'Q in Fig. ^(b) is then overlayed on 
the original region OQF. We see that any {m 2 i^ ym 2 i^) scatter plot will therefore exhibit 
two characteristic types of population density. For example, in the blue-hatched red area of 
Fig. ^(b) we expect the density of points to roughly double, since the folded distribution 
ON'P'Q is overlaid on top of the existing distribution OQF underneath. In Fig. || below, 
we shall mark such "double-density" regions with a blue cross-hatch in addition to the solid 
red shading. In contrast, region FQP'N' in Fig. ^(b) is a "single density" region, since the 
folded distribution happened to fall onto empty space, where originally there were no points 
to begin with. A single density region can also be obtained when portions of the original 
( m 'jln ,rn 'jif} sca tt er pl°t ar e not overlaid in the process of folding. In either case, we shall 
denote a single-density region by a solid (red) color-shading, but no cross-hatch. 

We are now ready to apply the intuition gained from Figs. ^ and [/] and identify the 
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Figure 8: The generic shape of the bivariate distribution versus m 2 -^^ for each of the four 

parameter space regions: (a) Region TZ\, (b) Region 7?.2, (c) Region IZ3 and (d) the off-shell case 
of Region H4. Each panel shows the typical shape (red-shaded) of the resulting (to^« \, m^,^) 
distribution, after the "folding" in Fig. ^. Blue-hatched (unhatched) areas correspond to double- 
density (single-density) regions. Each panel also shows the original location of the point P in the 
(m| ; , irfji ) plot, as well as the allowed positions of point P, following the color conventions of Figs. || 
and ^|. 

characteristic shapes of the ( m ^(/o)' m j«(?w)) distribution for each parameter space region 7£j. 
Our results are displayed in Fig. g, where we show the four characteristic shapes for each case: 
the on-shell cases of (a) Region 7£i, (b) Region IZ2, (c) Region 7^.3, and the off-shell case of (d) 
Re gion 7^.4. Each panel shows the typical shape (red-shaded) of the resulting (m^^ o y m 'ji^f l ^) 
distribution, after the "folding" in Fig. |?|. Blue-hatched (unhatched) areas correspond to 
double-density (single-density) regions. In addition, we show the original location of the 
point P in the (m^,m^) plot. The allowed positions of point P in each case are color- 
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shaded, following the color conventions of Figs. |2| and [6|. 

The nice feature of all the plots in Fig. || is that they are composed entirely of straight 
lines. This is a consequence of the fact that the original trapezoid in Fig. || is made up of 
straight lines, and then the "folding" of Fig. ^ does not curve the boundaries. Notice also 
the presence of internal kinematic boundaries, marking abrupt changes in the density of the 
distribution, e.g. O'N' in Fig. |(a), QHN' in Fig. |(b) and QF in Fig. |(c). 

It is clear from Fig. || that the shape of the ( m |;(/ )> scatter plot allows us to 

uniquely determine the parameter space region at hand. For example, the typical shape 
for Region TZi, exhibited in Fig. ||(a), consists of a right-angle triangular region OO'N 1 of 
double density and a right-angle trapezoidal region N'O'PF of single density. In this case, 
point P is directly observable, and its coordinates immediately yield the quantities n and p 
defined in ( |4.3| ) and (f4.5|). In addition, one can also measure the location / of point F along 
the irfjUfc} axis, given by eq. (4.4). This gives a total of three independent measurements: 
n, p and /, which should be ordered as n < p < f, in accordance with ( 4.15| ). Now we 
can clearly see the benefit of considering the two-dimensional ( m ^in o y m 'jUf l i0 distribution as 
opposed to the two individual one-dimensional distributions m 2 jin \ and rn~^ hi y Those one- 
dimensional distributions are obtained by projecting the ( m fln )i m "jl(hiO sca tter plot shown 
in Fig. |8](a) onto the two axes. It is easy to see from Fig. ||(a) that in this case the endpoint 
of the one-dimensional w |;n ) distribution will be given by c = n, while the endpoint of the 
one-dimensional m^u^A distribution will be given by d = /, and neither of those will reveal 
the quantity p. In contrast, p can be easily identified on the scatter plot, and provides an 
additional independent measurement. 

The case of Region 1Z%, which is shown in Fig. ||(c), is rather similar: the double-density 
region is still a triangle (OQF), while the single-density region is a quadrilateral (FQP'N 1 ). 
This time instead of point P we can clearly see its image P' , whose coordinates nevertheless 
still reveal the values of n and p. Point F is now hidden within the scatter plot, but may still 
be identifiable, since it corresponds to an abrupt change in density of points. Finally, now 
we have an additional measurement q of point Q, which is where the original line segment 
FP was folded in the reordering process of eqs. (L~l|,1.2). As a result, in region IZ3 we 
have a total of 4 measurements of kinematic endpoints, n, p, f and q, ordered as follows: 
p < q < f < n (see also eq. ( 4.17|) ). Later, when we project onto the two axes, the endpoint 
of the one-dimensional m j/(; ) distribution will be given by c = q, while the endpoint of 
the one-dimensional n^ufy^ distribution will be given by d = n. This now leaves out two 
additional potential measurements, p and /, which can be accessed on the two-dimensional 
scatter plot. 

In the third on-shell case of Region IZ2, shown in Fig. ||(b), the shape is more complex: 
the double density region is now a quadrilateral OQHN', while there are two disjoint single 
density triangular regions QP'H and N'HF. The point H appears on the intersection of the 
original FP boundary of the {m 2 ji n , m jif ) scatter plot in Fig. |6| and the (horizontal) image 
P'N' of the (vertical) PN boundary in Fig. ||. Once again, the coordinates of point P' reveal 
p and n, while points F and Q reveal / and q, correspondingly. In region 7Z2, therefore, 
there are 4 potential measurements, p, n, f and q, ordered as follows: p < q < n < f . When 
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the scatter plot of Fig. ||(b) is projected onto the axes, the endpoint of the one-dimensional 
m |i(/o) distribution will be given by c = q, while the endpoint of the one-dimensional m 2 ^^ 
distribution will be given by d = f. Once again, this leaves out two additional potential 
measurements, p and n, which can be extracted from the two-dimensional scatter plot. 

Let us now turn to the off-shell scenario of Region which is represented in Fig. ||(d). 
Because of the symmetry between the "near" and "far" leptons in the off-shell case, the folded 
region ONQ has an identical triangular shape as the underlying region OFQ, so that after 
the fold the two match perfectly and we obtain a single triangular region of double density, 
and no single-density areas. As can be seen from Fig. g(d), the off-shell scenario offers only 
one nontrivial endpoint measurement, which can be taken as /. The latter appears as the 
endpoint d in the one-dimensional n^^i) distribution, while the endpoint c of the other 
one-dimensional distribution, m 2 ( ^ o y is then simply given as c = //2, in agreement with 
eq. (^32|) . 



We have just seen that in the off-shell scenario of Fig. |l](b) the two-dimensional scatter 
plot (m^i(i j, m2 ji(hi)) does not yield any additional kinematic endpoint measurements. How- 
ever, it can still be helpful in discriminating a potential regional ambiguity which may arrise 
as follows. Notice that the triangular double-density shape of the scatter plot in the off-shell 
case of Fig. |8|(d) can in principle also be obtained in the on-shell cases of Fig. |(b) and 
Fig. |§|(c), provided that the image P' of point P ends up very close to point F. In terms 
of the (mji n ,m,ji ) scatter plot of Fig. ||, this situation corresponds to the on-shell cases of 
Regions IZ2 or TZ^, with point P lying very close to the yellow-shaded dot representing Re- 
gion 72-4. In spite of having the same shape of their boundary lines, the two scatter plots will 
be quite different, as they will exhibit a different point density. In particular, for all three 
on-shell cases, the pure phase space two-dimensional differential distribution ( |4.2| ) is given by 
the following (unit-normalized) formula 



d 2 T 



dm 2 jln dm 2 jlf n ^ m 2^ ( m 2j) fn-(f- p) m\ 



(for R BC < 1). (4.18) 



Notice that within the kinematically allowed region, the density is independent of m2 ji f - in 
the limit p — > 0, the expression ( 4.1S| ) becomes singular when m 2 i n — ► n. This singularity is 



regularized by the width of particle B and the branching fraction for the C — > B decay. In 
contrast, the corresponding density in the off-shell case is quite different, and in particular 
does not exhibit such singular features. 

We are now ready to revisit the duplication problem discussed in Section ||. We have 
just seen that the two-dimensional distribution of irijin ) versus fn 2 ^^ can help resolve the 
duplication in two very different ways. First, the shape of the kinematic boundary lines in 
the {fn 2 ui \i m2 ji{hi)) scatter plot uniquely identifies the region, as shown in Fig. |8|. Since the 
duplicate solutions that we found always appear in two different regions, this is in principle 
sufficient to eliminate the wrong solution. Secondly, the scatter plots offer the possibility of 
additional measurements, and at the very least a measurement of the quantity p. As can be 
seen from Table 0, the value of p is already different for each pair of duplicate spectra, and, 
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Figure 9: Scatter plots of f^m -\ versus m^,^ for the four study points from Table || exhibiting 
duplication: (a) point P31, (b) point P23, (c) point P32 and (d) point P23. Notice the quadratic scale 
used on both axes. The kinematic boundary lines are outlined with the corresponding color for each 
region, following the color coding conventions of Figs. || and ^. Each plot has 10,000 data points. We 
assume that all particles A, B, C and D are exactly on-shell. 



provided that it can be measured with sufficient precision, can also be used to remove the 
ambiguity. 

Our conclusions are confirmed by Fig. ^, which shows the ( m 'ji(i y rn ^i(f l i)^ sca ^ter plots 
for the four study points from Table |2| exhibiting duplication: (a) point P31, (b) point P23, (c) 
point P32 and (d) point P% 3 . The figure indeed shows that each pair of duplicate points has 
identical values for the endpoints of the separate one- dimensional invariant mass distributions 
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2 

m- 



• ( y and m 2 -^ hi y However, the shapes of the scatter plots are very different, and so are 
the values of the corresponding p endpoints. We therefore conclude that the duplication 
encountered in Section || ceases to be a problem, once we generalize the analysis to two- 
dimensional (bivariate) distributions as discussed here. 

In conclusion of this section, we point out that when the two-dimensional scatter plots 
like those in Fig. || are projected onto the axes to obtain the corresponding one-dimensional 
distributions of either m^in ) ov rn 2 },^, the latter often exhibit some peculiar features near 
their endpoints, which were classified as either "feet" or "drops" in ref. [30]. The origin of 
these features is now easy to understand in terms of the original two-dimensional scatter 
plot. For example, consider the scatter plots in Figs. ||(b) and ||(d). When projected onto 
the 

m ji(hi) ax i s ; both of them will exhibit a classic "drop" at the endpoint, which is 

simply due to the flat upper boundary P'N' in Fig. ||(c). Similarly, the projection of the 
scatter plot in Fig. |9|(c) onto the Tn 2 ^^ axis will exhibit a classic "foot" extending from n to 
/. The "foot" can be easily understood in terms of the generic shape of Fig. ||(b), where it 
arises from the projection of the single density area N'HF. 

5. Kinematic boundary lines for the vs. m 2 n distribution 

Following the logic of the previous section, we shall now proceed to analyze the two-dimensional 
distribution 

d 2 T 



dm 2 jU dm 2 



(5.1) 



whose generic shape OVUS is shown in Fig. 1C . The kinematic boundary lines of the scatter 



plot Q5.1D generally consist of four segments. The upper (SU) and lower (OV) curved bound- 
aries are parts of a hyperbolic curve OWS, while the left (OS) and right (UV) boundaries 
are straight lines. Therefore, in order to describe the shape of the (m^, m 2 u ) scatter plot, it 
is sufficient to provide the parametric equations for the upper and lower curved boundaries 
SW and OW, plus the location of the vertical line UV. In analogy with ( [4.9|) , we choose the 
variable on the horizontal axis, in this case m?,, as the line parameter describing the hyperbola 
OWS. Then the upper boundary line SUW is given by the parametric equation [33] 



m 2 M+) (ml) = l±^2± + l m l(l-R OD )(l-RAc) (5.2) 



I -Rod 



2 >| \ 

- 4:171 pR AC } 



while the lower boundary line OVW is given by [33] 

7 2 , 2n _ j+RcD m 2 1 
^u(-)^u) - 2 Rc^ + 2 



<u-M) = 1 -^^ + lm 2 D (l-R CD )(l-R AC ) (5.3) 



Rod 



Am%R A c 



(Tfb \ 
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Figure 10: The generic shape OVUS of the bivariate distribution ( p.l| ) in the (w ;i ,raj n ) plane 



The vertical straight line segment UV is in general located at mf, = a, where a is the value 
of the dilepton invariant mass endpoint (m™ ax ) 2 already introduced in Section 2T. As we 



saw in Section 2.1, the expression for a = (rrifi ax ) depends on whether we are dealing with 
the on-shell scenario of Fig. |l|(a) or the off-shell scenario of Fig. |l|(b). Therefore we shall now 
introduce separate notation for the endpoint a in each of these two cases. In the on-shell 
scenario of Fig. |l](a) we shall use a on to designate our previous eq. (|2.8|) 

a on = m 2 D R C D (1 - Rbc) (1 - Rab), (5-4) 



while in the off-shell scenario of Fig. |T|(b) we shall use a Q g for the previous result fl2.18| ) 

flofr = m 2 D R C d (1 - \/Rac) 2 - (5.5) 



From these two equations, it is not difficult to see that 



(5.6) 



as indicated in Fig. 10. As an aside, we mention that with the help of eqs. (|5.4j ) and 
it is easy to show that the equal sign in (|5.6|) is achieved when Rab = Rbc, i- e - when the 
on-shell spectrum happens to lie exactly on the border between regions Tl 2 and K 3 in Fig. g. 
In physical terms this means that the mass of particle B is equal to the geometric mean 
of tua and mc- 

tub = ^m,Am c => a on = a oS . (5.7) 

This represents another potential source of confusion in extracting the mass spectrum - 
the measurement of the dilepton invariant mass endpoint a alone tells us nothing about 
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whether the intermediate particle B is on-shell or off-shell, leading to two possible solutions 
[28]. Fortunately, with the inclusion of the additional measurements c, d and e, we did not 
encounter this type of duplication in the course of our analysis in Sec. ||. 

Returning now to our discussion of the (mf^ m 2 u ) scatter plot in Fig. [TIJ eq. ( |5.6| ) implies 
that in the on-shell case of Fig. 0(a), the data points do not fill up the whole region OWS, 
but only extend up to the vertical boundary UV. The region to the right of the UV line is 
kinematically inaccessible. On the other hand, in the off-shell case of Fig. |(b), the whole 
region OWS is filled up. The only exception to this rule is the very special on-shell case of 
( |5.7| ), when the UV line moves to the very tip W of the hyperbola, thus allowing the whole 
region OWS, as if this were an off-shell scenario. 

In analogy with our discussion in Section ||| of the kinematic boundaries in the two- 
dimensional (jrijin \-i m2 ji(hi)) distribution, we now identify several special points along the 
hyperbola OWS in Fig. [l0|. Point O is simply the origin (0,0) of the (rri 2 l ,rn 2 ll ) coordinate 
system. Point W is the tip of the hyperbola, where the upper branch m^^^m^) meets the 
lower branch m^^_-j(m^). By definition, the coordinate of point W is a a g, while its m 2 u 
coordinate is 

w = m ^(+)(«off) = raf H (_)(a ff) =m 2 D (l- Rcd\/Rac) U - y/ Rac) ■ (5-8) 

Point S is where the upper kinematic boundary line m^^^m^) intersects the coordinate 
axis. The rn 2 u coordinate of point S is therefore 

* = "*?«(+) (0) =m 2 D (l- Rod) (1 - Rac) ■ (5.9) 

Points U and V label the intersections of the vertical boundary UV with the upper and lower 
hyperbolic branches ( |5.2|) and (j5^3|), respectively. They share the same m| coordinate a on , 
while their coordinates are correspondingly given by 

u = m2 iiK+)( a on) (5-10) 



1 2 

2 D 



(1 + R CD )(1 - Rbc){1 - Rab) + (1 - Rcd)(1 - Rac + \Rbc ~ Rab\ 



V = m2 ill(-)( a on) (5.11) 



1 2 

2 D 



(1 + Rcd)(1 ~ Rbc)(1 - Rab) + (1 - r cd){1 - Rac - \Rbc - Rab\) 



Finally, there is one more special point on the upper branch SUW: it is the point T where 
minimization condition 

/ dm 2 ., 



m %u+)( m 'ii) h as a local maximum. The m| coordinate at of point T can be found from the 



L m+) 



dm l L* 



= (5.12) 

mf l= at 



and is given by 



at = m 2 D [RcD-y/RAD)n-y/RAD). (5-13) 
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Then, the m 2 - u coordinate t of point T is easily found by substituting (|5.13|) into (p^): 

t = m %(+)( a t) = m 2 D (l - yfR^Y . (5.14) 



We should point out that point T as we have defined it here, does not exist in all parameter 



space regions. To see this, let us calculate the slope of the upper branch m^^.Amfj) at point 



dm? 



\ _ Rcd ~ Rac ^ 15 x 



dm u I 2 „ r cd (1 - Rac) 



Since the denominator is always positive, the sign of the derivative is determined by the 
relative size of Rcd and Rac- When Rcd < Rac, the slope is negative, and T does not 
exist. In that case, the maximum value of over the whole scatter plot OVUS is obtained 
exactly at S, and is given by s in eq. ([5.9|). Comparing to the first lines in eqs. ( |2.9| ) and 
( |2.19 ), we see that this happens precisely for the cases of Nju = 1 and Nju = 5. In contrast, 



for the other four cases Nju = 2, 3, 4, 6, the slope at point S is positive and point T is well 
defined. However, this does not mean that point T would then necessarily belong to the 
scatter plot OVUS. In the off-shell case of Nju = 6, point T clearly belongs to the scatter 
plot, and the maximum value of m| M is given by t in eq. ( |5 .14 ), in agreement with the second 



line in eq. (2.19). However, in the remaining three on-shell cases Nju = 2,3,4 one has to be 



more careful. Since the scatter plot is then limited by the UV vertical boundary, point T 
will be included only if it lies to the left of the UV line, i.e. we must have 

a t < a on . (5.16) 
Using ( 5.13j ) and ( |5.4[ ), this condition can be equivalently rewritten as 

(R BC - RabRcd){Rab - Rbd) > 0. (5.17) 
Alternatively, the point T will fall outside the scatter plot, whenever 

a t > a on , (5.18) 

or equivalently, 

(R BC - RabRcd){Rab - Rbd) < 0. (5.19) 
We see that whether point T is included or not, depends on the sign of the expression 

(R BC - RabRcd)(Rab - Rbd). (5.20) 

Notice that the two factors entering this expression cannot be simultaneously negative: if 
that were the case, we would have 



Rbc ~ RabRcd < => m 2 B <^m 2 c 
Rab ~ Rbd < =>• m A m D < m? B 



m D < mc, ( 5 - 21 ) 
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which contradicts our basic assumption ( |2.6| ). Therefore, whenever one of the two factors in 
( l5~20l) is negative, the other is guaranteed to be positive. Of course, it is also possible that 



both factors in ( |5. 20 ) are positive to begin with. Altogether, this leads to three different 



possibilities, which are related to the Nju = 2,3,4 cases of eq. ( p.9| ). 



Nju = 2. In this case, the first factor in eq. ( 5.20| ) is negative, leading to the following 
logical chain 

Nju = 2 : R BC - RabRcd < =► Rab - Rbd > 

=> (Rbc ~ RabRcd)(Rab - Rbd) < => a t > a on , (5.22) 

placing point T outside the scatter plot. Then, the maximum value of rn 2 u is obtained 



at point U and is given by eq. ( 5.10 ). Since in this case Rbc < RabRcd < Rab, the 
absolute value sign in ( 5.10| ) can be resolved as \Rbc — Rab\ = Rab — Rbc and then 
eq. (|5.10| ) simplifies to 

u = m 2 D (l - R BC )(1 - RabRcd), (5.23) 
confirming the result on the second line of eq. Q2.9| ). 

Nju = 3. In this case, it is the second factor in eq. fT2Cp which is negative: 



Nju = 3 : Rab - Rbd < Rbc - RabRcd > 

(Rbc - RabRcd)(Rab - Rbd) < a t > a on . (5.24) 

Once again, point T is outside the scatter plot, and the maximum value of is 
obtained at point U and is given by (|5.10| ). This time, however, Rab < Rbd = 
RbcRcd < Rbc, and correspondingly, \R B c ~ Rab\ = Rbc ~ Rab- Then, eq. (|5.1C| ) 
simplifies to 

u = m D {l-R AB ){l-RBD), (5.25) 



agreeing with the third line of eq. 
Nju = 4. This is the case when both factors in eq. ( |5.20| ) are positive, leading to 

Rbc - RabRcd > 1 
N jU = 4: \ (R BC -RabRcd) (Rab -Rbd) > => a t < a on . 

Rab - Rbd > J 

(5.26) 

Point T now belongs to the scatter plot, and its coordinate t defined in ( 5.14| ) gives the 



maximum value of the distribution, in agreement with the fourth line of (|2.9| ). 



Note that Fig. |10| now allows us to understand geometrically the physical meaning of 
the lower threshold e = (^^jj^jl)) 2 introduced in Section ||. If we restrict ourselves only to 
points with > ^a on , i.e. to the right of the dashed line EE', the one-dimensional m| H 
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Figure 11: The same as Fig. ||, but for versus m^ u . 



distribution will exhibit a lower endpoint, whose value e is given by the m 2 u coordinate of 



point E in Fig. 10. In the on-shell case, e is given by 



e = mJ H( _ ) (a ou /2) = -m 2 D i (1 + R CD ){1 - Rab)(1 ~ Rbc) 



(5.27) 



+ 2(1 - R CD )(1 - R AC ) - (1 - RcdW(1 + Rab) 2 {1 + Rbc) 2 - IGRac }, 
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while in the off-shell case e is given by 



= m 2 jll( _ ) (a oS /2) = Kn 2 D (l - y^)^2R CD (l - y/W^) (5.28) 
+ (1 - Rcd) U + VRao -\]l + R A c + G^/Rac) I ■ 



It is not difficult to see that eqs. ( 5.27| ) and ( p.28 ) are identical to ( |2.15| ) and ( p. 22 ), corre- 
spondingly. 

The newly introduced quantities s, t, u and v can be directly observed experimentally 12 
on the scatter plot of Fig. Table || lists their square root values for our four duplicate 
study points P31, P23, -P32 and P^. As expected, the value of u is matched identically for each 
pair. However, the two other directly observable quantities s and v differ, and in principle 
can be used to resolve the duplication. This is illustrated in Fig. |ll], where we plot the two- 
dimensional distribution fl5.1[ ) for each duplicated example: (a) P31, (b) P23, (c) P32 and (d) 
P23 • Unlike Fig. ||, here the differences between the scatter plots for each duplicated pair are 
only quantitative, and may be difficult to observe in practice. The fact that the plots look 
similar is not very surprising, given our earlier discussion. Notice that duplication occurs 
only in regions with Nju = 2 or Njn = 3. In both cases, the shape of the (m 2 ;,m| w ) scatter 
plot is rather similar: the slope at point S is positive, and the upper boundary SU is cut 
off before it reaches the local maximum at point T. Furthermore, the duplication analysis 
ensures that the rightmost vertical boundary UV occurs in the same location Q on . 



6. Summary and outlook 

We now summarize the main results of the paper and discuss possible directions for future 
investigations. Our main results are as follows: 



In Section |2.2| we provided analytical inversion formulas which allow the immediate 
calculation of the mass spectrum rriA, itib, mc an d mo hi terms of a set of four 
measured invariant mass endpoints {a, c, d, e}. Our formulas are valid in all parameter 
space regions, since we do not use the endpoint b = m™ x , which is problematic in 



regions (3,1), (3,2) and (2,3), see eq. (2.17). 



Once the endpoint m^ x is eliminated from the discussion, we only need to consider 4 
different cases, IZi, i = 1, 2, 3, 4, as illustrated with the color-coded regions in Fig. ||. In 
contrast, previous studies which made use of the m 1 ^^ endpoint [26,27,30] were forced 
to consider all 11 different possibilities (Nju,Nji) shown in Fig. |2| 



12 



In principle, one can also measure indirectly the locations of points W and T, even when they fall outside 



the observable scatter plot. Since the analytical expressions (5.2) and (5.3) for the boundary lines are already 
known, one can fit them to the observable portions on the scatter plot, and then extrapolate the obtained 
analytical fit into the kinematically inaccessible region, thus obtaining the "would-be" positions of T and W . 
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• We investigated analytically the possibility of finding multiple solutions for the mass 
spectrum, even when a perfect experiment can measure the values for all five invariant 
mass endpoints {a, b, c, d, e} with zero error bars. Although we still had to consider sep- 
arately each of the four different cases IZi , we found that in most of the parameter space 
the spectrum is uniquely determined. Unfortunately, there is also a certain portion of 
parameter space, illustrated in Fig. [|, where one finds an exact duplication, i.e. two 
very different mass spectra yield identical values for all five measurements {a, b, c, d, e}. 
The situation is only going to get worse, once we take into account the inevitable exper- 
imental errors on the endpoint measurements, which can only proliferate the number of 
candidate solutions. Our results show that the conventional method of invariant mass 
endpoints may not be sufficient and one needs to look for new fresh ideas. 



The main goal of this paper is to advertise a new approach to the study of the usual 
invariant mass distributions. In particular, we point out that the multivariate invari- 
ant mass distributions contain a lot more useful information than the individual one- 
dimensional histograms, which are usually considered. As two illustrative examples, 
we discussed the two-dimensional {i r n^m \ > m ^t{hi) } distribution in Section || and the 
two-dimensional {rri^, m 2 u } distribution in Section [|. The former is always bounded 
by straight lines (see Fig. p|), while the latter is bounded by the hyperbola given by 



(|5.2|) and (|5.3|), and (in the on-shell case only) by the straight line UV in Fig. 10. 



The two-dimensional distributions exhibit two useful features. First, their shapes, 
i.e. the locations and orientations of their boundary lines, are characteristic of the 



corresponding parameter space region TZ{, as shown in Figs. |§| and 1C. This observa- 
tion can be used to identify the relevant parameter space region, and resolve potential 
ambiguities in the extraction of the mass spectrum. Second, the boundary lines exhibit 
a number of special points, whose coordinates can in principle be measured, providing 
additional experimental information about the mass spectrum. For example, in the 
■f m j£(Zo)' m j£(hi)} sca tter plots of Fig. || one may identify points F, P, N and Q, and 
correspondingly measure their n r >,jUf li \ coordinates, which (in the on-shell case) are given 
by 



(6.1) 
(6.2) 
(6.3) 

(6.4) 



/ = 


m D 


(1 


- R CD ) (1 - Rab) , 


p = 


m 2 D 


(1 


- R CD ) Rbc (1 - Rab) , 


n = 


m 2 D 


(1 


- R C d) (1 - Rbc) , 


Q = 


m 2 D 


(1 


1 - Rab 
~ Rcd) 2-R ab - 



Similarly, on the m 2 u } scatter plot in Fig. 10 one may identify the points S, U, 
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V, E, and (sometimes) T and W. Their coordinates are given by 
s = m 2 D (l- R CD ) (1 - Rac) , 

(1 + R C d)(1-Rbc)(1-Rab) 



m 2 D 



(6.5) 
(6.6) 



1 2 

2 D 



+ (1 - R C d)(1 - Rac + I^bc - Rab\) 



(1 + R CD )(1-Rbc)(1-Rab) 



+ (1 - /2ou)(l - Rac - \Rbc - Rab\) 



(6.7) 



1 + Rcd)(± - Rab)(1 ~ Rbc) + 2(1 - Rcd)(1 - Rac) (6.8) 
-(1 - RcDWi 1 + RabYO- + RbcY - 16i?AC , 



t = m 2 D [l- i/R 



IAD 



1 - iicDx/^) (l " V^Ac) 



(6.9) 
(6.10) 



The advantage of the new approach is apparent from eqs. ( 6. l - |6.ld| ). Including the 
dilepton invariant mass endpoint a, the set of potential invariant mass endpoint mea- 
surements has now expanded to 11: 



{a, f, P, n, q, s, u, v, e, t, w} 



(6.11) 



instead of the original five: 

{a, b, c, d, e}. (6.12) 

Of course, the endpoints in Q6.11| ) are not independent from each other, since they 
are all given in terms of only 4 input parameters fl2.7|) . Nevertheless, it is certainly 
preferable to have as many measurements as possible. The redundancy of information 
is a virtue, since it helps to improve the precision of the mass determination. 

The inversion formulas may simplify considerably, if we replace e, whose analytical ex- 
pression fl6.8p is rather complicated, with some of the other measurements in (6.1 - |6 . 1 0] ) . 
One such example is shown in Appendix [A|, where we start from the set {a, /, p, n}, 
and obtain a very simple result ( |A.ll| - [ATl4|) for the inversion. 

An important advantage of the two-dimensional approach is that one can readily resolve 
the ambiguity between the endpoints of the and the m 2 iL distributions. Indeed, 



"jln 



notice that the endpoints (6.1) and Q6.3|), are region- independent, and can be directly 
observed from the boundary lines. This removes the need to consider the different 
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parameter space regions IZi one by one. The possibility of distinguishing the jl n and 
jlf invariant mass endpoints from two-dimensional scatter plots was also suggested in 
Refs. [88] and [89], where the {mf^m^} distribution was used instead. 

• Another advantage of the two-dimensional representation of the data is that one can 
then perform a fit to the boundary lines of the scatter plot instead of a fit to the 
endpoints in the one-dimensional distributions. This improves the precision of the 
mass determination, as demonstrated in [89] for the SPSla SUSY benchmark example. 

In conclusion, we outline several directions for future investigations. 

★ Perhaps the most pressing question is whether and how well the method proposed here 
will survive the experimental complications of a full-blown analysis including detector 
simulation, backgrounds from Standard Model as well as SUSY combinatorics, the finite 
widths of the particles B, C and D, the varying population density of the scatter plots, 
etc. This is currently under study in the CMS SUSY working group and results will be 
presented in a separate publication. 

* In this paper we limited ourselves to the analysis of the boundary lines of two-dimensional 
distributions. However, the method can be easily generalized by including one more 
dimension and studying the boundary surface of the three-dimensional distribution 
( |1.9[) . A similar generalization was already shown to be beneficial in the case of spin 
measurements [69]. 

* One could also consider other choices of two-dimensional distributions, for example 

{ m h m jl}> { m % m2 ju} I 88 !' or { m «> m J(k>)}> { m U' m2 ji(hi)} I 89 ]- Those distributions 
also allow to discriminate between the "near" and "far" lepton endpoints, and will 
contribute even more data points to the set ( 6.11| ). 



-k One could also generalize the method to a longer decay chain, e.g. one which starts 
with a gluino [27]. 
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A. Appendix: Simple inversion formulas in regions TZ\, IZ2 and IZ3 

In Sec. U have saw that the shape analysis of a { m2 ji(i y m ?ji(hi)) scatter plot alone reveals the 
values of 

/ = KT) 2 = m 2 D (l - R CD )(1 ~ Rab), (A.l) 
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P = (m^lY = m 2 D R BC (l - Rcd)(1 - Rab) = fRsc, (A.2) 

n = {mfC? = m 2 D (l - R C d)(1 - Rue), (A.3) 

in each of the three on-shell regions IZi, IZ2 and IZ3. In addition, in regions IZ2 and 7^3 one 
also has a fourth measurement 

9 = ( m S) 2 = m » C 1 " ^ ^f^f • (A- 4 ) 
Given these four measurements, it is worth asking whether the spectrum of four masses 
rriA, tub, tuc and thd can be uniquely determined based on the (^jm )i m ^i(fii)) sca ^ ev pl°t 
alone. Unfortunately, this is not possible, since the four measurements /, p, n and q are not 
all independent, due to the constraint (4.8). Therefore, one more independent measurement 
is needed. 

Fortunately, the dilepton mass edge measurement is both robust and on-shell-region- 
independent. Thus adding 

a = (mD 2 = m 2 D R CD (l - Rbc)(1 ~ Rab), (A.5) 

we obtain a set of 4 measurements 

{a,f,p,n} = {(mj^f, (mg) 2 , {m^f, KT) 2 }, ( A -6) 

which can be easily inverted to obtain the spectrum: 

Rab = 1- — , (A.7) 
n 

Rbc = j, (A.8) 

Rod ={l + ^) 1 , (A.9) 

In terms of the actual masses we get 

2 an P f-, f-P 



2 an P 1 A io\ 

m * = (/- p )2 ' (A - 12) 



an f 

(F^pf 



m c = 77 — 7^2' (A.13) 



m D = (/_ p )2 ( 1 + — ) • ( A - 14 ) 



Notice the simplicity of these formulas in comparison to ( 2.23j -2.2C) and ( 2.2^ - p730| ). The 
simplicity is mostly due to the fact that we are not using the measurement fl2.15| ) whose 
analytical expression is rather complicated. 
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